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THE  CRREL  HOPK.INSON  BAR  APPARATUS 
Piyush  K.  Dutta  Dennis  FarreLl  and  John  Kaiafut 


INTRODUCTION 

The  safe  design  of  structures  at  low  temperatures  must  consider  the 
brittle  behavior  of  material  under  dynamic  loads.  In  cold  regions,  numer¬ 
ous  modes  of  dynamic  loading  on  structures  are  encountered:  wave  action, 
wind  action,  ice  movement  and  impact,  drilling  and  machinery  vibrations, 
vehicular  motion,  pile  driving,  blasting,  and  earthquakes  are  but  a  few 
examples.  The  hulls  of  ice-breaking  vessels  are  subjected  to  repeated 
impact  loading  by  ice.  Therefore,  sound  engineering  design  for  cold 
regions  must  consider  material  property  data  determined  at  low  temperature 
and  at  high  strain-rate  loading.  The  CRREL  Hopkinson  bar  apparatus  (HBA) 
has  been  designed  to  perform  such  high  strain-rate  loading  tests.  This 
report  describes  the  design  and  operational  procedures  of  this  apparatus 
and  the  technique  used  to  analyze  the  data.  In  this  technique,  deformation 
characteristics  of  materials  at  high  rates  (50  to  1000  strains  per  second) 
are  obtained  by  analyzing  the  stress  waves  through  a  test  specimen  sand¬ 
wiched  between  two  elastic  bars. 

To  illustrate  the  capability  of  the  testing  method  some  experimental 
results  are  given. 


BACKGROUND 

The  name  of  the  Hopkinson  bar  originates  from  the  work  of  the  British 
physicist  B.  Hopkinson  (1914),  who  carried  out  impact  tests  on  various 
materials.  In  his  tests  he  generated  compressive  stress  pulses  in  long 
bars  by  impacting  them  on  end  with  bullets.  The  compressive  pulse  is 
reflected  at  the  opposite  free  end  of  the  bar  as  a  tensile  stress  pulse, 
and  brittle  materials  such  as  rock  or  concrete  fracture  in  tension  under 
the  influence  of  this  tension  pulse.  Kolsky  (1949,  1953)  first  used  the 
idea  of  stress  wave  propagation  from  one  bar  to  another  through  a  “est 
material  sandwiched  between  the  two  bars  and  thus  gave  the  method  wider 


applications.  This  method,  now  called  the  split  Hopkinson  bar,  involves 
the  determination  of  dynamic  stresses,  strains,  or  displacements  occurring 
at  the  end  of  a  bar  by  various  sensors  monitored  some  distance  away.  The 
documentation  of  the  apparatus,  instrumentation,  and  analysis  is  quoted 
frequently  in  the  literature  (Llndholm  1964,  Zukas  et  al.  1982). 

In  the  conventional  configuration,  the  apparatus  is  used  mostly  for 
compression  tests.  However,  various  researchers  have  modified  the  appara¬ 
tus  to  conduct  tensile  tests  (Harding  et  al.  1960,  Hauser  1966,  Llndholm 
and  Yeakley  1968,  Christman  et  al.  1971,  Nicholas  1981,  Harding  1983,  Ross 
et  al.  1984).  Duffy  et  al-  (1971)  described  using  the  split  Hopkinson  bar 
to  achieve  high  strain  rate  in  torsional  loading.  The  primary  advantage  of 
the  torsional  mode  of  wave  propagation  is  its  nondlspersiveness.  In 
addition,  the  radial  inertia  effects  are  not  present.  Nicholas  (1975) 
extended  the  technique  to  conduct  a  dynamic  three-point  bend  test  (Charpy 
impact  test)  using  a  single  bar  instrumented  with  a  strain  gauge. 

Klepaczko  (1980)  has  suggested  the  use  of  the  Hopkinson  bar  method  to  study 
the  dynamic  fracture-initiation  properties  of  materials.  He  used  a 
wedge-loaded  specimen  in  a  conventional  split  Hopkinson  bar  arrangement  to 
obtain  dynamic  load-displacement  data. 

THEORY  OF  OPERATION 

A  schematic  of  the  CRREL  Hopkinson  bar  used  for  compressive  testing  is 
shown  in  Figure  1.  The  striker,  driven  by  compressed  air,  approaches  from 
the  left  and  impacts  the  incident  bar.  This  sets  up  a  compressive  stress 
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Figure  1.  Split  Hopkinson  pressure  bar  system  sche¬ 
matic. 
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Figure  2.  Lagranglan  x-t  dia¬ 
gram  of  stress  pulse  wave  pro¬ 
pagation  in  split  Hopkinson 
pressure  bar. 

wave  pulse  in  the  incident  bar.  The  amplitude  of  the  stress  pulse  depends 
on  the  impact  velocity,  and  the  duration  depends  on  the  length  and  material 
characteristics  of  the  striker  and  the  incident  bar.  The  compressive 
stress-wave  pulse  propagates  through  the  incident  bar  and  arrives  at  the 
incident  bar/specimen  interface.  The  details  of  the  wave  propagation  in 
the  bars  are  illustrated  with  the  Lagranglan  x-t  diagram  shown  in  Figure 
2.  At  the  incident  bar/specimen  interface,  part  of  the  incident  pulse  is 
reflected  because  of  the  mismatch  of  the  mechanical  impedances  of  bar  and 
specimen.  The  rest  is  transmitted  into  the  specimen,  propagates  through 
it,  and  reaches  the  spec iraen/transmitter  bar  interface.  At  this  second 
interface,  part  of  the  stress  pulse  is  also  reflected  back  into  the 
specimen  and  the  remainder  is  transmitted  to  the  transmitter  bar.  If  the 
specimen  is  short,  the  wave-transmitting  time  will  be  small  compared  with 
the  duration  of  the  loading  stress  pulse.  Thus,  many  wave  reflections  take 
place  within  the  specimen  during  a  single  pulse  reflection,  so  stress  and 
strain  along  the  specimen  can  be  assumed  to  be  uniform. 

Consideration  of  equilibrium  at  the  interface  between  specimen  and 
transmitter  bar  shows  that  the  force  in  the  specimen  and  in  the  transmitter 
bar  is  equal-  The  force  is  measured  by  measuring  the  strain  on  the  elastic 
transmitter  bar  using  strain  gauges.  At  the  location  of  the  strain  gauge 
this  force  will  have  a  time  shift  represented  by  the  travel  times  from  the 
interface. 

The  average  strain  Eg  in  the  specimen  is  calculated  from  the  dis¬ 
placements  at  the  end  of  the  specimen.  The  theory  is  well  documented 
(Zukas  et  al.  1982),  and  it  is  summarized  in  Appendix  A  for  both  short  and 
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long  specimens, 
shows: 


For  short  specimens  the  analysis  presented  In  Appendix  A 


Specimen 

strain 

es(t)  = 

r  e  (t)dt 

(1) 

spec imen 

strain  rate 

dE^(t)  _ 

•  -2c  ,  .  , 

s 

(2) 

dt 

specimen 

stress 

a  (t)  = 
s 

E^et(t) 

(3) 

where 

c  =  wave  propagation  velocity  In  the  Hopklnson  bars 
Lg  =  length  of  specimen 
ej.(t)  =  Instantaneous  reflected  strain 
Cj-Ct)  =  Instantaneous  transmitted  strain 
E^j  =  elastic  modulus  of  the  bar  material. 

Equation  1  shows  that  the  average  strain  In  the  specimen  can  be  com¬ 
puted  by  measuring  the  reflected  strain  In  the  incident  bar.  This  is  made 
possible  by  locating  the  strain  gauge  a  sufficient  distance  from  the  inter¬ 
face.  This  allows  the  Incident  stress  pulse  to  pass  through  this  strain 
gauge  station  completely  before  the  arrival  of  the  reflected  pulse  from  the 
interface.  Thus,  a  longer  pulse  requires  a  correspondingly  longer  incident 
bar.  A  longer  pulse  Is  desirable  because,  as  mentioned  earlier,  the  strain 
uniformity  within  the  specimen  under  repeated  reflection  can  be  ensured. 
Therefore  longer  bars,  within  practical  limits,  are  more  convenient  and 
provide  more  accurate  measurements  than  shorter  bars. 

Equation  2  shows  that  for  high  strain  rate  the  reflected  wave  Ej. 
should  be  large  or  the  specimen  short.  Note  also  that  a  larger  reflection 
happens  only  when  there  is  a  larger  mismatch  in  impedance*  between  the  bar 
and  the  specimen.  This  occurs  only  when  the  specimen  is  of  smaller  dia¬ 
meter  than  the  bar  or  the  bar  is  stlffer  and  of  higher  density  than  the 
specimen  material- 

Areas  of  cross  section  (A),  modulus  of  elasticity  (E) ,  and  density  ( o) 
are  the  parameters  that  control  the  Impedance.  For  example,  the  impedance 
of  body  A  matches  the  Impedance  of  body  B  if  =  Ig,  i.e. 

.  -0.5  .  ^0.5 

*  Mechanical  Impedance  of  an  Interface  between  two  longitudinal  bodies 
determines  the  wave  propagation  characteristics  through  that  interface. 


If  >  Ig ,  a  compressive  stress  wave  reflects  back,  as  a  tensile  wave 
from  the  Interface,  whereas  If  <  Ig  a  compressive  stress  wave 
reflects  back  as  a  compressive  wave  and  vice  versa. 

Testing  of  crystalline  or  gram '.a r  materials  usually  requires  that 
specimens  be  5  to  10  times  (preferably  10  to  20  times)  the  maximum  grain  or 
crack  size.  Thus  larger-dlameter  specimens  are  preferable.  Again  the  dla- 
me ter-to-length  ratio  (2:4)  makes  the  specimens  in  these  categories  of 
materials  longer  than  the  metal  specimens  tested  In  the  Hopkinson  bars. 
Another  disadvantage  of  granular  material  is  the  comparatively  lower  wave 
velocity  within  the  specimen.  Thus,  for  a  relatively  short  incident  wave¬ 
length,  the  number  of  internal  reflections  within  the  specimen  Is  not  large 
enough  to  make  it  reach  a  state  of  uniform  stress.  Under  such  conditions 
longer  incident  wavelength,  and  therefore  longer  incident  bars,  are  desir¬ 
able. 

With  the  electronic  recording  and  computational  devices  available,  it 
is  possible  to  record  digitally  all  three  waves  (l.e.  incident,  reflected, 
and  transmitted  wave),  superimpose  them  one  over  the  other,  and  analyze  the 
displacements  of  the  two  interfaces  of  the  specimen.  The  average  stress 
and  strain  of  the  longer  samples  can  be  computed  as  shown  below  (refer  to 
Appendix  A  for  derivation  of  these  formulas): 


average  specimen  strain  e  (t)  =  — - —  J  [o. (t)  -  a  (t)  -  a  (t)  Idt  (4) 

s  pc  Lj  i  r  t 


average  specimen 
strain  rate 


average  specimen  stress  \  <^^(1)  crj.(t)] 


It  is  obvious  that  the  above  approach,  to  study  the  dynamic  properties 
of  materials,  is  indirect.  The  strain  rates  that  occur  during  the  test  are 
primarily  a  function  of  the  test  material's  own  constitutive  behavior.  It 
is  also  important  to  note  that  the  stress,  strain,  and  strain  rate  are 
average  values  and  are  calculated  on  the  basis  of  the  uniaxial  stress-state 
ass  ump  1 1 o  n  s . 


DESIGN  DETAILS 

The  Hopkinson  bar  apparatus  (Fig.  3)  Installed  at  CRREL  Includes  three 
major  subsystems:  the  stress  generating  system,  the  specimen  cooling 
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the  bars  must  be  properly  aligned.  Proper  alignment  is  also  Important  to 
ensure  that  minimum  bending  waves  are  generated.  Otherwise  these  waves  are 
sensed  by  the  strain  gauges  and  distort  the  true  uniaxial  compression/ 
tension  signals.  Each  bar  is  guided  through  three  support  blocks.  The 
support  blocks  and  the  launch  cylinder  are  mounted  on  a  long  rigid  table. 
Locations  of  the  blocks  and  the  impactor  are  adjustable  to  ensure  adequate 
alignment  during  test.  To  reduce  dispersion  of  the  stress  waves  through 
the  support  members,  the  support  blocks  were  specially  designed  with 
low-friction  ball  bearings  (Fig.  4).  Adjustments  are  provided  on  each 
individual  block  so  that  the  striker,  incident  bar,  specimen,  and  trans¬ 
mitter  bar  are  all  in  alignment. 

Hie  striker,  the  incident  bar,  and  the  transmitter  bar  are  made  of  the 
same  material  --  303  stainless  steel.  The  stainless  steel  bars  were  chosen 
for  their  noncorrosive  characteristics.  On  impact  the  stress  developed  in 
the  bar  should  be  held  below  its  yield  stress  [67,000  Ib/in.  "  (462  MPa)]  . 

To  limit  the  stress  level  to  50%  of  the  yield  stress,  it  was  calculated 
that  the  striker  velocity  should  not  exceed  38  ft/s  (11.6  m/s).  Calculat- 
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Figure  5-  HPB  apparatus  striker. 


ing  the  kinetics  of  the  striker  motion  shows  that  at  100  Ib/in. ^  (689  kPa) 
pressure  a  2-in.  (51-mm)  stroke  is  sufficient  to  reach  this  velocity. 

The  striker  is  12  in.  (305  mm)  long  and  has  the  same  diameter  as  the 
incident  bar  and  the  transmitter  bar  [1.5  in.  (38  mm)].  Both  the  incident 
bar  and  the  transmitter  bar  are  8  ft  (2438  mm)  long  and  have  ends  machined 
normal  to  the  bar  axis.  Tl-ie  striker  (Fig.  5)  end  is  rounded  with  a  spheri¬ 
cal  radius  of  12  in.  (305  mm)  to  minimize  any  misalignment  between  the 
striker  and  the  incident  bar. 

■b  -3 

To  cause  the  striker  to  impact,  an  8.5  ft  (0.24  m  )  compressed-air 
gas  reservoir  is  first  charged  to  the  required  pressure.  This  has  been 
calibrated  to  give  the  desired  stress  level  in  the  incident  bar.  Figure  6a 
shows  the  air  pressure  and  the  corresponding  incident  stress  level  data.  A 
quick-operating  solenoid  valve  is  then  activated  to  allow  the  compressed 
air  to  drive  the  striker,  which  impacts  the  incident  bar.  Figure  6b  shows 
the  repeatability  of  the  stress  wave  forms  with  six  separate  impacts. 

The  striker  itself  is  housed  in  a  launch  cylinder.  The  striker  bear¬ 
ing  surface  and  the  cylinder  inside  surface  form  a  close  fit,  with  a  0.003- 
in.  -(0  .08-mm-)diaraeter  clearance  between  the  striker  and  the  cylinder. 

This  minimizes  air  leakage  and  provides  guidance  for  the  striker.  A  0.25- 
In. -( 6 .4-mm-)d iameter  rod  connected  to  the  rear  of  the  striker  bar  pro¬ 
trudes  from  the  launch  cylinder  and  provides  a  convenient  means  to  locate  the 
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a.  At  various  air-pressure  levels  used  tor 
driving  the  striker. 
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b.  Repeatability  of  the  pulses  under  sIk 
separate  striker  impacts  at  20  Ib/in.^ 


Figure  6.  Stress  wave  pulse  forms. 


striker  inside  the  launcli  cylinder  and  adjust  the  stroke  length-  The  air 
in  front  of  the  striker  is  vented  to  atmosphere  through  two  0 . 75-in. -( 19- 
mm-)diaraeter  holes  in  the  front  part  of  the  launch  cylinder.  In  its 
current  configuration  no  attempt  has  been  made  to  measiire  the  velocity  of 
the  striker  at  impact.  This  is  considered  a  redundant  measurement  becaus 
the  c-l(;ctric  strain  gauging  is  available  to  measure  accurately  the  stress 
levi.'l  in  the  incident  bar  directly.  The  inside  surface  of  the  launcti 


cylinder  is  hoiiod  to  re^.ace  friction  and  wear.  To  reduce  friction  further, 
the  striker  bearings  are  inade  of  Teflon. 


The  specimen  is  located  between  the  incident  bar  and  the  transmitter 
bar.  A  s:rnll  amount  of  preloading  is  necessary  to  hold  the  cylindrical 
specimen  in  position  between  the  two  bars.  This  is  achieved  with  two 
rubber  bands  tensioning  the  two  bars  to  close  on  the  specimen.  In  addi¬ 
tion,  two  Teflon  collars  mate  the  ends  of  the  specimen  to  both  bars.  Each 
collar  is  slit  lengthwise  to  allow  a  modest  flexibility  in  diametrical 
expansion  when  required  for  mounting  the  specimen. 

Tlie  striker  together  with  the  incident  bar,  transmitter  bar,  and  the 
shock  absorbing  device  constitute  a  set  to  test  material  at  stresses  less 
than  33,500  lb/ in.  (23i  MPa). 

The  specimen  materials  tested  to  date  are  brass.  Teflon,  and  poly- 
crystalllne  ice.  The  next  section  describes  the  results  of  testing  an  ice 
Siimplo . 

Tlie  ice  specimens  were  prepared  in  molds  by  using  snow  particles  of 
sizes  0.07  to  0.03  in.  (1.78  to  0.833  mm).  The  mold  used  to  prepare  the 
samples  is  shown  in  Figure  7.  It  was  first  lubricated  and  tempered  at  10°F 
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(-12'’C).  Then,  before  adding  and  tamping  the  prepared  snow  in  the  three 
ctiambers,  the  system  was  tempered  again  at  O^F  (-18°C).  A  distilled,  de¬ 
sired  water  supply  was  gravity-fed  through  the  bottom  end  caps,  and  the 
slight  overflow  was  drained  through  holes  in  the  top  end  caps  to  remove  any 
trapped  air.  With  the  water  still  slowly  circulating,  the  mold  was 
returned  to  lO^F  (-12°C)  fo’’  freezing.  A  Styrofoam  insulation  cap  was 
applied  to  the  top  of  the  mold  to  promote  initial  ice  growth  from  the 
bottom  of  the  mold  and  additional  air  escape  at  the  top. 

After  the  specimens  were  frozen,  they  were  tempered  several  minutes  at 
40'’F  (4.4°C)  to  free  the  ice  from  the  mold.  After  they  were  ejected,  each 
specimen  was  wrapped  in  cellophane  and  stored  at  10°F  (-12°C). 

Specimen  Cooling  System 

To  keep  the  test  specimens  cold  during  the  impact  test,  the  specimen 
and  the  mating  ends  of  both  bars  are  enclosed  in  two  copper-colled  enclo¬ 
sures  through  which  cooled  nitrogen  gas  is  circulated.  The  nitrogen  gas  is 
cooled  by  circulating  it  through  a  liquid  nitrogen  bath.  The  details  of 
the  cooling  system  are  shown  in  Figure  8.  The  coils  themselves  are 
enclosed  in  a  Styrofoam  container.  The  specimen  temperature  is  monitored 
with  a  thermocouple. 

Temperature  is  controlled  by  the  level  of  immersion  of  the  cooling 
coil  in  the  liquid  nitrogen  and  the  flow  rate  of  the  nitrogen  gas  through 
the  coils.  With  this  arrangement  a  stable  temperature  as  low  as  -90°F 
(-68°C)  can  be  achieved.  In  the  current  design  configuration  the  N2  gas 
flow  is  branched  to  two  specimen-cooling  coils.  Continuous  operation  of 
the  coils  also  chills  the  ends  of  the  Hopkinson  bars,  permitting  rapid 
restali lizatlon  of  specimen  temperatures. 

Stress  Measurement  Systetn 

In  performing  the  experiments  it  is  necessary  to  record  the  complete 
profiles  of  the  incident  stress  wave  pulse  and  reflected  stress  wave  pulse 
in  the  incident  bar  and  the  transmitted  wave  pulse  in  the  transmitter  bar. 
For  this  purpose  foil  strain  gauges  were  bonded  on  both  the  incident  and 
transmitter  bars  at  the  midpoint  of  each  length.  The  positions  of  the 
gauges  and  the  lengths  of  the  bars  were  so  selected  that  the  stress  wave 
signals  could  be  recorded  for  their  entire  duration  without  interruptions 
caused  by  wave  reflections  from  the  ends.  Two  active  strain  gauges  were 
mounted  on  each  bar,  one  diametrically  opposite  the  other  to  cancel  bend¬ 
ing  wave  strain,  if  any  is  produced  because  of  misalignment. 
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Figure  8.  Coi:)ling  system:  spec  Imeti'cooling  coil  circulating  cooled  N2  gas 

The  schematic  for  the  eloctroulc  circuitry  is  shown  in  Figure  9.  The 
strain  gauge  outputs  after  prearapiif ication  are  fed  to  a  Model  4562  plug-in 
unit  oi  a  N'icolet  4(i94A  Digital  Oscilloscope.  Channel  A  is  used  for 
recording  the  incident  and  reflected  waveform  and  channel  B  for  the  trans- 
:riitted  wave.  The  waveforms  are  sampled  at  0.5-iis  intervals  and  digitized 
with  12-bit  resolution.  Tlie  scope  has  the  capability  to  expand  the  wave 
forms  after  recordi:ig.  It  can  store  and  recall  data  from  a  floppy  disk. 

Tlie  scope  is  compatible  with  the  Hewlett  Packard  Model  7470A  X-Y  plotter 
tor  plotting  the  displayed  waveforms  on  paper. 

An  on-line  data  acquisition  and  analysis  system  has  been  assembled 
using  a  PCFyhHb  interface  board  (mannf acttired  by  Capital  Equipment  Corp.) 
in  an  IBM  PC/AT  comiiuter.  Suitable  software  has  been  developed  to  trans- 
mi  t  I  h.-  St.  )r.,-d  or  recorded  data  from  the  scope  to  the  computer  for 
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Figure  9.  Stress  pulse  waveforms  recording  electronic 
circuit  schematic. 


In  initiating  the  test  the  sweeps  were  triggered  by  electrical  contact 
between  the  striker  and  the  incident  bar.  Figure  9  includes  the  schematic 
of  the  triggering  circuit. 

The  system  was  calibrated  dynamically  by  stress  wave  recording  induced 
in  both  bars  butted  together  without  a  specimen.  Figure  10  shows  these  two 
waveforms.  Note  that  on  the  first  trace  the  reflected  wave  from  the  butted 
surface  is  only  4%  of  the  total  energy  transmitted  to  the  second  bar.  In 
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Figure  10.  Incident,  reflected,  and 
transmitted  stress  pulse  wave  from  a 
single  impact  on  the  bars  belted  to 
each  other. 
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Figure  12.  Incident,  reflected, 
and  transmitted  stress  pulse  waves 
reconstituted  on  the  computer  screen. 


the  oscilloscope  screen.  Figure  12  shows  the  hardcopy  printout  after  the 
data  are  reconstructed  by  the  computer  program.  The  HOPP  program  then 
continues  to  analyze  these  waveforms  using  eq  4 ,  5,  and  6  to  produce  the 
time-variant  data  on  specimen  average  strain,  specimen  stress,  and  strain 
rate.  The  stress-strain  data  is  displayed  on  the  CRT  screen  and  can  also 
be  printed  out  on  a  printer.  The  plot  of  the  strain  rate  vs  strain  is  also 
displayed  on  the  CRT  screen  and  again  can  be  printed  out.  Figure  13  shows 
the  plot  of  the  stress-strain  curve  of  the  specimen  during  passage  of  the 
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stress  wave  and  Figure  14  shows  strain  rate  vs  strain.  If  desired,  a 
point-by-polnt  data  analysis  report  at  specified  time  intervals  in  tabulat¬ 
ed  format  can  also  be  printed  out. 
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Figure  13.  Dynamic  stress- 
strain  relationship  of  an  ice 
sample  computed  from  the  stress 
pulse  waveforms. 


Figure  14.  Strain  rate  vs. 
strain  of  an  ice  sample  during 
passage  of  the  stress  wave. 
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GENERAL  COMMENTS 


The  Hopklnson  bar  apparatus  with  Its  associated  data  recording  and 
analysis  equipment  is  a  convenient  means  to  perform  a  quick  evaluation  of 
materials  under  high  strain  rate  loading.  However,  it  must  be  recognized 
that  the  approach  is  an  indirect  one,  and  the  strain  rates  that  occur  dur¬ 
ing  the  test  are  primarily  a  function  of  the  material's  own  constitutive 
property.  Moreover,  two  Important  assumptions  are  made  in  the  analysis 
that  influence  the  results:  1)  nondispers i vi ty  of  the  elastic  waves  in  the 
bar  and  2)  the  absence  of  radial  inertia. 

For  these  reasons  the  split  Hopkinson  bar  has  been  the  subject  of  an 
extensive  analytical  study  to  evaluate  the  validity  of  the  results  obtained 
from  its  tests.  Davies  (1948)  was  the  first  to  exaiaine  critically  the 
effects  of  wave  dispersion  in  the  Hopkinson  bar,  he  concluded  that  the 
dispersion  effects  can  disguise  the  stress  change  that  happens  within  any 
1-us  Interval.  Inertial  ef'jcts  were  studied  by  Davies  and  Hunter  (1963), 
who  recommended  Chat  to  minimize  error  due  to  inertial  effects  and  inter¬ 
facial  friction  the  optimum  geometry  for  the  specimens  should  be  a/h  =  1.15 
where  a  is  the  radius  and  h  is  the  thickness  of  the  specimen.  Interfacial 
friction  problems  were  studied  in  great  depth  by  Bertholf  and  Karnes 
(1975),  whose  results  show  that  lubricated  Interfaces  have  the  minimum 
error. 


Another  important  consideration  in  performing  the  Hopkinson  bar  test 
on  specimens  at  other  than  ambient  temperature  is  the  influence  of  Che 
temperature  gradient  along  Che  bar.  Malvern  (1984),  as  reported  by  Zukas 
et  al.  (1982),  has  considered  this  problem.  His  analysis  is  based  on  the 
assumption  that  the  modulus  E  is  linearly  dependent  on  temperature. 

Assuming  that  the  temperature  gradient  is  linear,  a  more 
straightforward  derivation  of  the  correction  factor  is  possible.  The 
relation  between  the  stress  near  the  sample  (og)  and  the  stress  measured 
at  the  gauge  station  (og)  can  be  obtained  by  considering  the  general 
relationship  (refer  to  Zukas  1982) 
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where  E,  is  the  modulus  of  elasticity  of  the  bar  material  at  the  gauge 
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station  at  ambient  temperature  T^,  and  Eg  is  the  modulus  of  elasticity 
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of  Che  bar  material  near  Che  specimen  at  test  temperature  Tg .  For  the 
CRREL  Hopkinson  bar  made  of  AISI  303  stainless  steel,  Eg  at  0“F  (-18°C)  = 
29.4x10'"  Ib/in.’  (202.6  GPa)  and  Eg  =  28.3x10®  (195.0  GPa)  at  -300°F 
(-L84"C).  Using  these  values,  the  error  due  to  a  0  to  -300°F  (-18  to 
-184"C)  temperature  differencial  will  be  approximately  2%. 

In  interpreting  the  Hopkinson  test  bar  data  it  is  also  important  Co 
take  into  consideration  the  stress  wave  reflections,  stress  nonuniformity, 
aiid  Large  variation  of  strain,  especially  during  Che  initial  and  final  por¬ 
tions  of  the  test.  The  stress-strain  data  are  thus  valid  only  after  some 
iegree  of  uniformity  in  stress  and  strain  is  achieved.  Figure  14  shows 
that  Che  strain  rate  changed  rapidly  till  Che  strain  was  established  at 
about  36  scrains/s.  Once  strain  rate  uniformity  is  established,  the 
dynamic  stress-strain  cliaracter isCic  data  are  more  valid. 

A  series  of  tests  tias  been  conducted  using  ice  samples  and  a  few  other 
engineering  materials.  These  data  are  Che  subject  of  another  report  now  in 
preparation.  However,  it  has  been  generally  observed  that  at  the  low 
stress  Level  in  a  semibrittle  material  such  as  ice,  Che  stress-strain 
curve  is  smoother  chan  Chat  obtained  with  higher  stress  level  loading 
(Fig.  15).  At  the  low  stress  level  only  a  few  cracks  were  induced  in  the 
sample,  whereas  at  a  high  stress  level  impact  the  specimen  shattered  (Fig. 
16).  Multiple  reflections  of  the  stress  waves  within  the  sample  during  the 
fracturing  may  have  contributed  to  the  sawtooth  characteristic  of  the 
stress-strain  curve  for  high-level  stress  loading. 

The  heart  of  the  system,  which  gives  an  indirect  measurement  of  the 
dynamic  constitutive  relationship,  is  the  correct  recording  of  the  waves  in 
the  pressure  bars.  The  three  stress  waveforms  --  oj^(t),  o^(t),  and 
^(.(t)  --  form  tile  very  basis  of  this  analysis.  Since  "i{(t)  is  recorded 
earlier  and  both  '’t(t)  are  recorded  later  than  they  occur  at 

the  interface,  the  recorded  stress  waveforms  miist  be  appropriately  shifted 
in  time  bef(3re  they  are  combined  as  in  eq  4,  5,  and  6.  Tlie  use  of  a  digi¬ 
tal  oscilloscope  and  the  data  transfer  and  analysis  techniques  now 
dev.'loped  1 1 1  ow  the  Incident  bar/specimen  interface  and  transmitter  bar/ 
specimen  interface  displacement  to  be  determined  accurately  --  thus  the 
iverage  strain  in  the  specimen  can  be  computed  by  taking  the  arithmetic 
mean  of  Lhe*two  interface  displacements.  The  deviation  from  the  average 
d  i  spl  acemeiU  is  estiraaticd  as  half  the  difference  between  the  two. 
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This  procedure  has  aoc  been  used  much  in  the  past.  Most  of  the  data 
in  the  literature  were  obtained  before  automated  data  recording  and  analy¬ 
sis.  The  usual  procedure  has  been  to  assume  that  tl.e  two  interface 
displacements  are  equal,  neglecting  strain  variation  along  the  length  of 
the  specimen.  For  short  specimens,  where  strain  uniformity  could  be 
achieved  rapidly  and  the  stress  pulse  is  of  long  duration,  this  assumption 
may  not  lead  to  any  great  error,  but  the  current  design  allows  a  better 
averaging  technique  for  comparatively  longer  specimens. 
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Figure  15  Typical  dynamic  stress-strain  data  from  tests  of  polycrystal- 
llne  ice  samples. 
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APPENDIX  A.  THEORY  OF  THE  APPARATUS 


Elementary  plane  wave  propagation  theory  (Zukas  et  al.  1982,  Malvern 
1984) . 

In  the  one-dimensional  propagation  of  a  stress  wave  through  a  solid 
bar  the  following  assumptions  are  made; 

1)  A  plane  cross  section  remains  plane  as  the  wave  propagates  along 
the  longitudinal  axis. 

2)  Stress  distribution  is  uniform  across  the  section  of  the  bar. 

3)  Radial  inertia  can  be  neglected. 
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V  1 

1 

A  |A' 
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Particle  velocity. V  =  ^  =  £C 


Figure  Al .  Elementary  uniaxial 
plane  wave  propagation. 


In  Figure  Al  consider  an  element  AB  having  a  length  1  on  this  bar.  A 
stress  pulse  of  amplitude  a  (and  strain  e  =  n/E)  traveling  with  a  wave 
velocity  c  arrives  at  plane  A  at  the  Instant  t  =  0.  It  arrives  at  plane  B 
after  a  time  1/c.  During  this  time  the  element  AB  is  under  compression  n, 
therefore  plane  A  moves  to  A'  such  that 
AA'  =  deformation  =  c  1  . 

The  velocity  of  A  (called  the  particle  velocity  V)  is  then  given  by: 


V  = 


AA' 


1/c 


1/c 


=  EC 


However,  it  is  well  known  that  in  an  elastic  medium  of  density, 
propagation  of  a  longitudinal  wave,  is  given  by 
.0.5 


(Al) 
D,  the 


c  =  (E/n) 


(A2) 


o  r 
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Since  '^  =  E  r  from  eq  A1  and  A2 


->  V 

t  =  c"o*--  =  ncV 
c 


(A3) 


the  displacement  U  of  the  particle  at  any  time  t  is  given  by 

t  1  t 

U  =  '"  V(t)dt  =  —  (  o(t)dt  (A4) 

0  0 


Split  Hopkinson  bar  wave  mechanics 

Consider  x  and  y  as  the  ends  of  two  bars  of  the  same  diameter  and 
material.  Any  wave  in  which  the  particle  moves  to  the  right  is  assumed 
to  move  in  a  positive  (+)  direction,  and  one  that  moves  to  the  left  assumes 
a  negative  (-)  direction.  Thus,  a  compression  wave  that  propagates  to  the 
right  has  the  particle  velocity  in  the  positive  direction,  and  a  tensile 
wave  that  moves  to  the  left  has  a  positive  particle  velocity. 

A  compressive  stress  wave  pulse  arriving  at  the  x-y  boundary  is 
partially  transmitted  and  also  partially  reflected. 

If  n|(t)  =  stress  wave  incident  to  the  sample  sandwiched  between  the 
two  bar  edges  x  and  y,  then  displacement  of  the  x  face,  due  to  the 

incident  stress  wave  pulse  o^Ct),  is  given  by 


U. 

1 


aj^(t)dt 


and  displacement  Uj.  of  the  x  face  due  to  the  reflected  stress  wave  pulse 
'''r(t)  is 

,  t 

U  =  —  ^  a  (t)dt 

r  pc  „  r 


As  the  reflected  wave  travels  in  the  direction  opposite  to  the  incident 
wave,  the  total  displacement  of  the  x  face  is  given  by 
,  t 

U.  -  U  =  ~  r  la.(t)  -  a  (t)]dt  .  (A5) 

1  r  pc  ^  1  r 

Similarly,  displacement  of  the  y  face  due  to  the  transmitted  stress  wave 
pulse  n(-(t)  is  given  by 
,  t 

U  =  '■  a  (t)dt  (A6) 


Therefore,  total  detonnatlon  U-p  ot  the  sample  due  to  passage  of  the 

stress  wave  pulse  through  it  is 

L’  =  U  .  -  U  -  U 
T  i  r  t 


"i'--  ' 

T  3C  g 


^iCt)  -  o^(t)  -  a^(t)]dt  . 


Note  that  to  obtain  U(-  from  eq  A6  ,  nj(t)  should  be  measured  at  the  same 
instant  that  np(t)  and  aj-(t)  are  measured.  If  ^t  =  time  required  for 
the  wave  to  propagate  from  face  x  to  y  and  Oj-Ct)  is  measured  at  time  T 
after  start  of  the  sweep,  then  a^-Ct)  should  be  measured  at  time  T  +  At. 
If  Lg  =  length  of  the  sample,  then  the  instantaneous  sample  strain 


CgCt)  is  given  by 


£  (t)  =  -- 
s  L 

s 


C  [ a. (t)  -  a  (t)  -  a  (t) ]dt 
1  r  t 


Therefore  the  sample  strain  rate  (tg)  can  be  obtained  by  differentiating 
the  above  integral  expression  Cg  with  respect  to  time,  i.e.  deg(t)/dt 
is  given  by 

r  =  — —  [a  (t)  -  0  (t)  -  n  (t))  .  (A9) 

s  ocL  1  r  c 

s 

Sample  stress  at  interface  x  is  given  by  the  algebraic  sum  of  the  incident 
stress  np(t)  and  the  reflected  stress  a^Ct);  i.e. 


the  stress  at  interface  x,  <7j^(t)  =  a^(t)  +  o^(t) 
the  stress  at  interface  y,  n^(t)  =  a^(t) 


(A9.1) 

(A9.2) 


The  average  sample  stress  of  ag(t)  is  given  by  the  mean  of  the  stress  at 
the  X  and  y  Interfaces,  i.e. 


a^(t)  =  ^  [Oj^Ct)  +  a^(t)  +  0^.(1)] 


(AlO) 


By  using  eq  A8 ,  A9 ,  and  AlO  the  complete  stress-strain  wave  and  the  strain 
rate  loading  in  the  HPB  set-up  can  be  determined. 
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Short  speci men 

If  the  length  of  the  sample  is  short,  then  the  stress  at  interface  x 
can  be  assumed  to  be  same  as  the  stress  at  interface  y,  i.e.  from  eq  A9.1 
and  A9 . 2 


= 

+  a^(t) 

or 

1 

f-T 

II 

(All) 

Subs t i tu t ing 

the  above  in  eq  A8 , 

t  (t)  = 
s 

-1 

t 

f  2a  (t)dt  . 

0  ^ 

(A12) 

DcL 

s 

Again, 

a^(t) 

=  E  G^(t)  =  C-OG^(t) 

and 

Cs(t) 

=  ^^-r  e(t)dt  ; 

therefore 

=  f  e  (t)dt 

^  /-I  ^ 

s  0 

(A13) 

and  strain  rate 

dt 

• 

r  = 

S 

-2c  ,  , 

-  e  ( t )  . 

L  r 

(A14) 

Again  substituting  eq  All  In  eq  AlO,  the  sample  stress  for  a  short 
sample  is  given  by 

i  [2a|.(t)]  =  fT^(t)  ,  (A15) 

t  herefore 


(A16) 


APPENDIX  B.  INSTRUMENTATION  OF  THE  HPB  APPARATUS 


Strain  was  measured  at  the  midpoint  on  each  half  of  the  split 
Hopkinson  bar.  Two  0.125-in.  (3.18-mm)  long,  350-ohm  strain  gauges,  com¬ 
pensated  for  stainless  steel,  were  attached  at  each  location.  One  was 
aligned  diametrically  opposite  the  other  for  bending  moment  compensation. 
The  gauges.  Type  CEA-09-125  UR-350,  were  bonded  with  AE-10  epoxy  and  cured 
at  room  temperature  according  to  accepted  procedures.  The  gauge  dimensions 
(0.125  X  0.60  in.)  were  chosen  so  strain  averaging  would  occur  over  only  a 
small  physical  area  of  the  bar. 

The  gauges  were  connected  into  a  bridge  configuration  with  the  three- 
wire  technique.  Each  gauge  cable  was  individually  shielded  and  the  gauge 
was  wired  into  opposing  sides  of  a  two-active-arm,  four-arm  bridge.  This 
configuration  offers  twice  the  voltage  output  of  a  two-gauge,  one-active- 
arm,  four-arm  bridge  with  less  nonlinearity  and  only  a  minimal  loss  of 
bending  compensation.  Micro  Measurement  M  Coat  D  was  used  over  the  com¬ 
pleted  gauge  Installation  for  moisture  protection.  Each  gauge  was 
connected  to  the  bridge  completion  and  balance  circuit  by  15  ft  (4.6  m) 
of  Belden  8723  cable  where  the  four-arm  bridge  was  completed  with  350-ohm 
fixed  resistors. 

Bridge  balance  was  achieved  with  the  500-kohm  potentiometer,  as  shown 
in  Figure  9  of  the  text.  The  balance  circuit  was  designed  as  described  by 
Stein  (1967).  Bridge  excitation  was  set  at  15V  and  was  provided  by  the 
Hewlett  Packard  HP  6218B  power  supply.  Each  bridge  had  its  own  independent 
power  source. 

Both  amplifiers  were  set  on  a  gain  of  ten,  which  provided  a  corre¬ 
sponding  bandwidth  of  100  kHz  with  120  dB  of  common  mode  rejection  at  their 
differential  Inputs.  The  output  of  each  amplifier  was  single-ended  and 
connected  into  the  positive  input  terminal  of  each  scope  channel.  Data  was 
recorded  on  the  +100  mV  settings  and  digitized  at  the  rate  of  0.5  us  per 
point,  giving  a  record  length  of  about  4  ms.  Two  channels  on  the  scope 
were  used:  channel  I  recorded  incident  and  reflected  strain  and  channel  2 
recorded  the  transmitted  strain.  The  oscilloscope  trace  was  triggered  when 
the  hammer  contacted  the  first  bar  of  the  split  Hopkinson  bar.  This  con¬ 
tact  applied  a  3V  battery  voltage  directly  to  the  external  trigger  input  of 
the  scope. 
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Figure  Bl-  Bridge  con¬ 
figuration  for  bending 
wave  effort  analysis. 


Figure  B2.  Bridge 
configuration  for 
shunt  calibration. 


The  bridge  configuration  used  for  recording  the  strains  in  the  bars  is 
shown  in  Figure  Bl.  In  this  model  the  strain  generated  by  bending  will 
have  an  insignificant  effect  on  the  longitudinal  strains  computed  from 
^out  • 

For  example,  In  the  bending  configuration  the  change  of  strain  ^R  will 
be  positive  in  one  gauge  and  negative  in  the  opposite  gauge,  so  that 

_ _ E(R+AR) 

^out  ”  =  2R  -  AR  "  2R  +  AR  ’ 


which  on  simplification  gives 


E 

out 


E 


4R^  -  Ar2 


(Bl) 


Since  AR  has  a  very  small  value  compared  to  R,  Eq^,.  for  bending  is 
negligibly  small. 

Calibration  was  performed  with  the  usual  gauge  shunt  technique.  When 
shunted  with  the  349,650-ohm  resistor  shown  In  Figure  B2  ,  each  active  gauge 
gave  the  equivalent  of  250  microstrain  at  the  oscilloscope  input,  as  shown 
be  low : 


R 


AD 


^  (350)(3496  50) 
R  +  Rg  350  +  349650 


349,650  ohm 


^i"ad 

R  +  r"^ 


340650  E 
350  +  349650 


0.49975  E 


E 

2  ■ 
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Therefore  Eq  =  -  e^  =  0.00025  E  or  E^/E  =  0.00025.  But 

,  _  gauge  factor  x  strain 

Eo/e  =  -2  , 

therefore 

0.00025  X  2 

strain  = - - - 

gauge  factor 

Assuming  the  gauge  factor  is  2,  strain  =  0.00025  =  250  microstrai n. 

External  noise  pickup  was  kept  to  an  acceptable  Level  with  careful 
eanlpraent  layout.  The  metal  table  supporting  the  split  Hopkinson  bar  was 
tied  to  a  handy  electrical  conduit  and  this  became  the  electrical  common 
point.  As  mentioned  earlier,  the  gauge  cables  were  individually  shielded 
and  tied  at  the  far  end  to  the  electrical  common  as  were  all  enclosures  and 
shields.  Bridge  configuration  and  excitation  were  chosen  to  give  a  high 
electrical  output,  and  the  dc  amplifier  provided  120  dB  of  common  mode 
rejection  from  dc  to  60  Hz. 

Tlie  scope  was  capable  of  outputting  the  digitized  data  directly  to  a 
plotter  although  most  data  was  stored  on  floppy  disks  and  later  transmitted 
to  a  computer  for  further  analysis. 


Table  B1 .  Instrumentation  components. 
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Description _ _ _ Manufacturer 


Model  No. 


4  Strain  gauge 

4  Calibration  resistors 

4  Bridge  completing 

resistors 

1  Bridge  completion 

and  balance  unit 

2  Power  supply 

1  Graphics  plotter 

1  Interface  board 

1  Digital  oscilloscope 

with  dual  amplifiers 
and  dual  disk  drives 
1  External  scope  trigger 

(3V  battery) 

Gauge  epoxy 
Gauge  cable 


Micro  Measurements 
Micro  Measurements 

Micro  Measurements 


CEA-09-125  UR-350 
W-349650-02 

S-350-01 


In-house 

Hewlett  Packard  HP  6218B 

Hewlett  Packard  HP  7470A 

Capital  Equip.  Co.  PC0488 

Nicolet  Instrument  Co.  4099A 

4562 

XF44 

Union  Carbide  Corp.  W-356 

Micro  Measurements  AE-IO 

Belden  8723 


